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The Raman response of the metallic glass NieyZra^ is measured as a function of polarization 
and temperature and analyzed theoretically. Unexpectedly, the intensity in the range up to 
300 cm~^ increases upon cooling, which is counterintuitive when the response originates from vi- 
brations alone as in insulators. The increase finds a natural explanation if the conduction electrons 
are assumed to scatter on localized vibrations with a scattering probability proportional to the 
Dcbye- Waller factor. None of our assumptions is material specific, and the results are expected to 
be relevant for disordered systems in general. 
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I. INTRODUCTION 

The intriguing properties of disordered or amorphous 
materials have been studied in great detail^ and play an 
important role in many contemporary studies and appli- 
cations. From the viewpoint of microscopic properties, 
one of the hallmarks is an excess density of states (DOS) 
in the vibrational spectrum as well as extra contributions 
to the low-temperature heat capacity and conduction and 
to the ultrasound attenuation^. The extra vibrational 
DOS which is superposed on the Dcbyo-like cj^ spectrum 
of crystals is usually referred to as the boson peak. In 
amorphous insulators a random distribution of the elas- 
tic constants is appropriate to capture the experimental 
results^. However, local modes, believed to be relevant 
for both disordered crystals and amorphous materials, 
can be the origin of the boson peak^. In metallic Ni- 
B, which exists in crystalline and amorphous form, extra 
vibrational states were observed in a neutron scattering 
study at both low and high energies''. In fact, the vi- 
brations in disordered systems can be phonons with con- 
siderable energies as well as low-lying excitations in the 
meV range usually present when one or a few atoms have 
more space and vibrate in voids, described as a poten- 
tial well. The latter excitations can be quite localized. 
The two lowest levels of these local modes are frequently 
called two level systems (TLS) and play an important 
role in the low-temperature thermodynamic behavior. 

Raman scattering experiments were instrumental in 
studying the vibrational properties of amorphous materi- 
als since the complete spectrum can be observed because 
of the absence of translational invariance and the related 
collapse of the selection rules. In an early calculation the 
corresponding Raman response x"j,(a;,r) was found to 
be proportional to the vibrational DOS g(uj) as^ 



where n(cj, T) = [(^'^/^bT _ ^-i jg ^-^^ Bose-Einstein oc- 
cupation mrnibcr at temperature T, S'i.s(cj) is the mea- 
sured Raman intensity (structure factor), and w = Wj— 
is the energy transferred to the system by incoming 
and scattered polarized photons i and s, respectively. 
The coupling factors c^^s, originally assumed to be con- 
stant, turned out to have a rather complicated energy 
dependence^'^. The factor w"^ is separated out of the 
coupling factor as the most probable energy dependence 
dnc to the matrix elements of the harmonic oscillators. 
Schmid and Schirmacher showed recently for insulators 
that g{ijj) is not simply the DOS. Rather, x"^(w,T) is 
a polarization dependent weighted sum of the longitudi- 
nal and the transversal susceptibilities in systems with 
fluctuating elasto-optical coupling and force constants.^ 
A similarly rigorous treatment for metals is still missing. 

In this paper we analyze Raman scattering results of 
the metallic glass NigyZraa and show that there are three 
relevant types of excitations which can be disentangled. 
There is a temperature-independent Drude-like response 
due to the electron-hole excitations as suggested earlier 
for dirty metals^. On top of the electronic spectrum, 
there is an additional response which, at 300 K, is strictly 
linear up to 20 cm^^ . In the range 20 to 300 cm^-' , there 
is a wide maximum which loses about 30% of its spectral 
weight upon increasing the temperature, just opposite 
what is expected from the occupation number of bosonic 
excitations. To our knowledge only the Debye- Waller 
factor may account for such behavior by reducing the 
scattering probability of electrons on vibrations. There- 
fore, substantial vibrational density of states is required 
at relatively low energies to produce a strong enough tem- 
perature dependence. We focus on this contribution. 

II. SAMPLES AND EXPERIMENT 



The amorphous NigyZraa sample was prepared by 
melt-spinning 99.9% pure Zr and Ni metals in an Ar at- 
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mosphere after electron-beam melting in a quartz tube. 
The non-crystalline state was checked by x-ray diffrac- 
tion. The amorphous character is also reflected in the 
slight (1%) increase of the resistivity p{T) upon cooling 
from room temperature to 4K, which is just opposite 
what is observed in microcrystalline samples^. A mean 
free path I of approximately 20 A or some 5 10 in- 
teratomic distances a is estimated from the resistivity, 
p = 150 /ill cm, and the density of states at the Fermi 
level N{EfY^. 

The Raman experiments were performed with stan- 
dard equipment. For excitation a solid state laser emit- 
ting at 532 nm (Klastech Scherzo-Denicafc-532-300) was 
used. The spectra are corrected for the sensitivity of 
the instrument. We always show the response function 
Rx'ls{^,T) cx S'.^s/jl -I- n{uj,T)}. R absorbs all exper- 
imental factors and s is the imaginary part of 
the dimensionless Raman response function. 



III. RESULTS AND DISCUSSION 
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In Fig. 1, the temperature dependence of the Ra- 
man spectra of Nig7Zr33 with parallel (VV) and per- 
pendicular (VH) light polarizations is shown. There 
is a slightly asymmetric maximum centered at about 
140 cm~^ (lmeV = 8cm~^ ) which is superposed on 
a linear continuum extending to several 100 meV 
[Figs. 1 (a) and 1(b)]. While the continuum is temper- 
ature independent the intensity in the peak maximum 
increases by approximately 30% upon cooling [Figs. 1 (c) 
and 1(d)]. 

The continuum has a strong resemblance to the be- 
havior of dirty metals for which the theory was worked 
out years ago*. In that work different scattering channels 
classified by crystal harmonics L were considered. The 
electron density fluctuations that appear in the channel 
L — are screened and cannot contribute to the spec- 
trum. However, the channels L — 2n > are relevant 
implying that the electrons are being redistributed be- 
tween different parts of the Fermi surface. Only those 
particles contribute to the cross section which are scat- 
tered out of a given channel. The scattering mechanism 
is similar to that in multi- valley semiconductors^^. In 
NigyZras the continuum can be explained essentially in 
this way with a relaxation rate in the range of 0.5 eV, 
consistent with the resistivity of 150 pft cm. The posi- 
tion of the low-energy peak would correspond to a much 
longer relaxation time (smaller rate). In addition, the 
spectral shape deviates from the expected one*. Hence, 
the peak has to be interpreted differently. 

Following the reasoning for insulators, one can argue 
that the peak resembles that of vibrational properties''^. 
Schmid and Schirmacher' make a prediction for the ra- 
tio of the responses in the VH and VV configurations, 
p{ijj, T) = Xy/z/xyy J derived for the case of locally vary- 
ing elastic constants. In our experiment, p{uj,T) in- 
creases with energy as shown in Fig. 1 (e) rather than 



FIG. 1: (Color online) Raman response Rx'{sjj,T) of 
Ni67Zr33 for perpendicular (VH) and parallel (VV) light po- 
larizations, (a) and (b) show on an extended energy range 
that the spectra consist of a polarization independent lin- 
ear continuum and superposed peaks in the range below 
300 cm~^ . (c), (d) After subtraction of the linear part, the 
low-energy peaks have nearly identical shapes for all tempera- 
tures [c.f. (e)], but the overall intensity in VH is only approx- 
imately 30% of that in VV. (e) The ratio p — Xvh/x'vv is 
temperature independent and increases slightly with energy. 



decreases as predicted. Whether or not the discrepancy 
can be traced back to the metallicity of NigrZrsa cannot 
be answered in the present state of knowledge. In ad- 
dition, in the case of a vibrational origin the intensity 
of the peak should increase rather than decrease with T 
since the occupation number of vibrations with energy uj 
is proportional to n{uj,T). In crystals with well defined 
electronic levels having widths in the range of fc^T, res- 
onance effects could reverse the temperature dependence 
expected for vibrations by overcompensating the effect of 
the Bose factor^^'^^. However, in an amorphous material 
with electronic relaxation rates in the range of eV, as di- 
rectly derived from the Raman experiment (see above), 
the resonances are too wide to be relevant. In fact, the 
electronic part is independent of the temperature. We 
have also verified experimentally that both the high- and 
low-energy contributions are independent of the energy of 
the exciting photons. Finally, the electrons could scat- 
ter from spin fluctuations which gain influence toward 
lower temperature in a paramagnetic solid. However, for 
Ni concentration below 70 % the paramagnetism is very 
weak^" as corroborated by our own magnetization mea- 
surements for the sample used here. Moreover, Batalla et 
al.^^ concluded from their studies of the low-temperature 
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specific heat that the coupling to spin fluctuations in Ni- 
Zr compounds is negligible in the composition range of 
interest here. Since, to the best of our knowledge, studies 
of the specific heat or of the vibrational spectra by other 
methods do not exist in the temperature range covered 
here further insight into the interpretation of the Raman 
results from a comparison to other work is not possible 
at present. 

It is instructive to divide the low-energy response 
[Figs. 1 (c) and 1(d)] by w as shown for the VH case 
in Fig. 2. In the VH spectrum at 315 K, which was mea- 
sured with high resolution, a constant part on top of the 
Drude response emerges below 20cm~^ as highlighted in 
the inset. In the spirit of Eq. (1) the constant part would 
correspond to a Debye-like variation of the DOS g{u!) 
if the coupling factors Ci^s were constant. However, in 
a large number of studies on insulators and semiconduc- 
tors, different energy dependences for Ci^s are found^'^. 
Further clarification of this question must be postponed 
to another publication. Above 20cm~^ , additional scat- 
tering sets in abruptly which, for the various reasons just 
discussed, cannot be interpreted in the same way as the 
boson peak in insulators. The increase toward low T in- 
dicates that other types of excitations must contribute to 
the inelastic relaxation of the electrons. 

In the following we study the electron relaxation on 
static defects such as impurities and on vibrations at en- 
ergies in the range of and below a maximal vibration 
frequency cj ~ on the order of the Debye energy. The 
thermal occupation number can affect the Drude peak 
and lead to an increase of the scattering amplitude in the 
spectrum with temperature. Given the observed temper- 
ature dependence this process cannot be the dominating 

one. In the range 20 300 cni^^ , an additional peak 

shows up on top of the Drude background. Electron-hole 
pairs are created which are coupled to local vibrations. 
If low-lying vibrations at w ^ wo are created the ampli- 
tude u of the vibrations can be large. We assume that 
u can be approximated by the result for the harmonic 
oscillator, \u{uj)\^ = {h/Mu}){l/2 + n{uj,T)}, where M 
is the relevant mass. Since the electron has a momen- 
tum close to the Fermi momentum, kp — 7r/a, the situa- 
tion is comparable to that of x-ray diffraction on crystals, 
and the scattering probability increases at lower temper- 
atures according to the Debye- Waller factor, e~^^^"'"^', 
with 2W{uj,T) = |k-up. In other words, the local os- 
cillator is a weaker scatterer at high temperature when 
it is blurred out because of its large amplitude and a 
strong one at low temperatures when it is almost point- 
like. Since n{uj,T) also appears in the exponent, the 
Debye- Waller factor could dominate over the occupation 
number at low energies as observed in the experiments 
(Fig. 1). On the other hand, static impurities and other 
imperfections of the system may be very rigid; thus their 
Debye- Waller factor is small. This is in accordance with 
the temperature-independent linear continuum. 

The Debye-Waller factor for a single vibration 
with energy can be estimated as 2Wj = 
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FIG. 2: (Color online) "Boson" peak in NigrZraa [same tem- 
peratures as in Fig. 1 (c-e)]. The inset shows the spectrum at 
315 K, which was measured down to 6cm~'^ , on a logarith- 
mic scale. The slight increase toward a; — >■ is a remainder 
of the laser line. At the other temperatures the spectra were 
measured only for cj > 15cm~^ . In the plot y^' /u a Drude- 
like spectrum appears as a constant (lower dashed line in the 
inset). The measured low-energy spectrum saturates above 
the Drude response (upper dashed line), indicating scatter- 
ing which could originate from vibrations with a Debye-like 
spectrum (see text). 

{a/ujj){\/2 + n{ujj,T)}, with a = 2P{h/M). M is the 
mass of the vibrating atom, and k ^ kp \{t: /a) where 
A is a number of order unity. Thus, a can be estimated as 
a ~ A\^LOp{ra/M) {ojp has the same units as ujj). With 
the magnitude of the Fermi energy Ep = hujp ~ 5 eV, a 
will be on the order 2A^ (cm~^). 

Now, the creation of a vibration at ujj by the electron- 
hole continuum in the presence of static impurities and 
disorder is considered*. The electron- hole continuum 
arising in channel L is described by the density operator 
Pi (9) = X]fc lL{k)a^ak+q and the corresponding correla- 
tion function is 7r° (w) = —2Np\^]^\'^ (1 — iujTl)~^, where 
Np is the density of states at the Fermi level, 7^ is the 
amplitude of the photon-electron- hole coupling and t£~^ 
is the relaxation rate. In the following we approximate 
t£ with To as they are not very different. The Hamilto- 
nian of the vibration is ujjb'^bj, where ojj, 6^, and bj 
are the energy and the creation and annihilation opera- 
tors of the phonon, respectively. The vibrational Green's 
function is given by the usual propagator multiplied by 
an exponential term due to the temperature dependent 
scattering probability, 

!?,(..) = -^e--("^-^-) (2) 
OJ -uj. 

The electron continuum is coupled to the vibrations by a 
local interaction, similar to that in small-polaron theory, 
as J2j 9L,jPLiRi) {b^ + bj), where gLj is the coupling 
and PL(.Ri) is the electron density operator at the posi- 
tion Ri of the local oscillator. The actual coupling is very 
model dependent but in order to demonstrate the anoma- 
lous temperature dependence we restrict ourselves to the 
simple Hamiltonian above. After averaging over the po- 
sitions of the static impurities tt^ becomes independent 
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of the position, and the renormahzed correlation function 



is 7ri(w) 



1 



The strength of the dimensionless couphng can be es- 
timated as NfJ2j91,jDj{u}) ^ NpNlVl'^/^j, where N 
is the number of the active vibrational centers per vol- 
ume, ujj is their typical energy (10 < ujj < 100 K) and 
V is the coupling potential. For a rough estimate typi- 
cal values are borrowed from the corresponding estimates 
for TLSsi^: \VNf\ < 0.2 and N « ^(e)wj, where N{e) 
is the energy density of the centers, e.g., N{e) « 5 x 
10~^ cm~^ = 5ppmK~^. Thus the order of magni- 
tude at the relevant energy window of that dimensionless 
coupling is NpNlVl'^/Cjj ~ N{e)\VNF\'^ /Np < 5 x lO'^. 
Assuming weak coupling the summation over the pertur- 
bation series is not required. Keeping first order only, 
ttlH « nliuj) l + Y.j9ljDj{uj)TTl{uj) , and the Ra- 
man response for the localized modes alone reads 

X gioc(w){l-f n(c.,T)}e-2^("'^) (3) 

and is plotted in Fig. 3. Equation (3) is valid in 
the limit lotq <C 1 and for uj > 0, where vibra- 
tions are created, and gioc(w){l + n(a;, r)}e"2'^("''^) = 

— TT^^/m jZ3j(w)|. Here, (7ioc(w) is proportional 

to the effective density of states of the local oscillators, 
which is weighted by the occupation number and the 
Debye- Waller factor. For uj > luq, gioc{(^) falls off, and at 
larger energies only the Drude peak survives. In the tem- 
perature dependence of the low-lying vibrational contri- 
bution, the Debye- Waller factor dominates over the tem- 
perature dependence of the occupation number n(w,T), 
and the low-energy part is suppressed as it is observed ex- 
perimentally. The residual scattering below 20cm~^ on 
top of the Drude response [Fig. 2 inset] will be studied 
in more detail in a future presentation. 

For the Debye-Waller factor to become effective, a 
substantial density of states at energies typically below 
100 cm~^ is required. The usual Debye-like spectrum 
with the DOS starting out as oj^ is inappropriate. TLSs, 
on the other hand, have excitations with extremely small 
energies. Since the temperatures of interest are relatively 
high we can mimic the potential as a single well with a 
harmonic spectrum and an occupation n ^ 1. In order 
to produce spectra we assume a constant .gioc(w) as the 
simplest example. We emphasize that there may be a 
substantial energy dependence, the discussion of which, 
however, is beyond the scope of this publication. Here, 
5ioc(w) is cut off smoothly at low and high energies by hy- 
perbolic tangent functions. As the low-energy vibrations 
are suppressed by the Debye-Waller factor, the results 
are insensitive to the low-energy cut off. In addition to 
the TLSs, other types of excitation may exist at very low 
energies. The actual calculation follows these steps: tq 
is fitted to describe the Drude peak at higher energies, 
and uiq is determined by the width of the vibrational 
peak. The effective density of localized states g\oc is con- 
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FIG. 3: (Color online) Response from localized oscillators 
Xjoc/^- The low-energy cut-off due to the Debye-Waller fac- 
tor is moving downward with decreasing temperature. At 
room temperature the cut-off can be clearly resolved in the 
experiments (inset of Fig. 2). 



stant with lower- and higher-energy soft cut offs at 40 
and 160 cm"^ with 20 and 90cm~^ widths, respectively, 
and / dujgioc = 1. Also, a = 80cm~^ is relatively large. 
However, the localized modes vibrate in voids with a very 
flat potential at equilibrium, and the harmonic approxi- 
mation underestimates the occupation number substan- 
tially. In addition, the momentum transfer k may be 
larger than tt/o, and the parameter A can easily be be- 
tween 2 and 4. 

Finally, the relative insensitivity of the spectral shapes 
at low energy to the light polarizations [see Figs. 1 (c)- 
1(e)] requires a comment since there is a discrepancy with 
the results of insulators^. In amorphous metals, due to 
the relatively short mean free path of the electrons, the 
intensities depend only on the light-electron bare ver- 
tex functions that contain matrix elements of the form 
{F\pes\v) {iy\pei\I) with /, F, v the initial, final and in- 
termediate states, respectively, p the momentum opera- 
tor, and e.i^s the photon polarizations. In the tight bind- 
ing approximation the intra-atomic matrix elements are 
the most important ones and can be transformed into 
dipole matrix elements. Hence, the selection rules for 
intra-atomic orbitals apply and can explain the weaker 
intensity of the VH component (crossed polarizations).^^ 
For completeness, we should mention that there is an 
additional peak at a few eV which is attributed to the 
multiple electron- hole excitations^''. 



IV. CONCLUSIONS 

In conclusion, we have shown experimentally and theo- 
retically that the low-energy part of the Raman spectra in 
the metallic glass NigyZras originates from the superpo- 
sition of electronic relaxation processes due to scattering 
from impurities and from localized vibrations. The tem- 
perature dependence of the spectral weight indicated that 
local modes contribute substantially to the cross section. 
The increase with decreasing temperature can in fact be 
described in a qualitative way assuming that the scat- 
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tcring of conduction electrons is dominated by the tem- 
perature dependence of the Debye- Waller factor rather 
than by the occupation number of the vibrations alone. 
For a quantitative comparison of experiment and theory, 
a concrete assumption for the effective density of states 
would be needed. As pointed out earlier, this would be 
very model dependent but would not influence the main 
conclusion concerning the influence of the Debye- Waller 
factor. The result may be applicable to a large variety 
of systems such as doped semiconductors, metallic func- 
tionalized carbon nanotubes, and polymers. 
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